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The 'elliptic cylindrical Kadomtsev-Petviashvili equation' 

( H T + 6HH^ + J^cCC + 7TT 1 ? 9\ H ~ 1 o ?/ ? ?\ ) H o kH uu — 0, 
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is derived for surface gravity waves with nearly-elliptic front, generalising the cylindrical KP equa- 
tion for nearly-concentric waves. We discuss transformations between the derived equation and 
two existing versions of the KP equation, for nearly-plane and nearly-concentric waves. The trans- 
formations are used to construct important classes of exact solutions of the derived equation and 
corresponding approximate solutions for surface waves. 

PACS numbers: 47.35.Fg, 47.35.Bb, 92.10.Hm 



I. INTRODUCTION 

There exist two classical versions of the Kadomtsev- 
Petviashvili (KP) equation pQ associated with the sur- 
face wave problems for an incompressible fluid, described 
by the full set of Euler equations with free surface and 
rigid horizontal bottom boundary conditions (see [2HI] 
and references therein): 



The original KP equation [1 

(U T + 6UU£ + Uztt)^ + 3a 2 U YY = (2) 
and the cylindrical KP (cKP) equation [2] 



W T + 6WW X + W xxx 



W 
2t 



3a' 



Wvv = (3) 



p(u t + UU X + VUy + wu z ) = —p x , 

p(v t + UV X + Wy + wv z ) = -p y , 

p(w t + uw x + vwy + ww z ) = -p z - pg, 

U X + Vy + W Z = 0, 

P\z=h(x,y,t) ~ Pa 

(1 + h 2 y )h xx + (1 + h 2 x )h yy - 2h x h y h xy 

" ^ (1 + ^ + ^)3/2 • 

™\z=h(x,y,t) =h t + uh x + v/fy, 



(1) 



Here, (u, v, w) are the three components of the velocity 
vector in Cartesian coordinates (x,y,z), t is the time, p 
the pressure (p a is the constant atmospheric pressure at 
the surface, and T is the coefficient of the surface tension), 
p is constant density, g is the gravitational acceleration, 
z = is the bottom, and z = h(x, y, t) is the free surface. 
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describe the evolution of nearly-plane and nearly- 
concentric waves, respectively. For surface gravity waves 
with no or weak surface tension <j 2 = 1 (KP-II), while 
the case a 2 = — 1 (KP-I) can be obtained when surface 
tension effects are large (see, for example, [5 and ref- 
erences therein). Transformations mapping solutions of 
the cKP and KP equations were independently found in 
[4j [6] . They have been used to construct some special 
solutions of the cKP equation in [7]. Indeed, the map 



W(t, X ,V)^U(t,£,Y) :=W[t,£ 



Y 2 



Y 



12cr 2 r' r 



transforms any solution of the cKP equation |3| into a 
solution of the KP equation (J2|. Conversely, the map 



U(t,Z,Y)^W(t, X ,V) :=U[t, X 



tV 2 
12a 2 ' 



-V 



transforms any solution of the KP equation Q into a 
solution of the cKP equation 

In this paper we derive a generalisation of the cKP 
equation (|3| for surface gravity waves, which can be writ- 
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ten in the form 



H r + 6HH C + tf ccc + 



2(r 2 



v 2 a 2 



12a 2 (r 2 -a 2 ) 



He 



3a 2 



— 0, 



where a 2 = ±1, describing waves with nearly-elliptic 
front. This elliptic cylindrical KP equation (ecKP) is 
derived from the full set of Euler equations for an incom- 
pressible fluid and free surface and rigid bottom bound- 
ary conditions ([!]), written in the elliptic cylindrical coor- 
dinate system. The linear long wave equation, written in 
these coordinates, does not allow for exact solutions de- 
scribing waves with elliptic front. However, there exists 
an asymptotic reduction to the necessary equation, and 
we show that this allows one to derive a generalisation of 
the cKP equation. 

We chose to derive the ecKP equation from the Euler 
equations rather than using the velocity potential formu- 
lation, since this opens the way to the study of internal 
and surface waves on a current for a fluid with arbitrary 
stratification, as well as accounting for the effects of a 
variable background and Earth's rotation (see [8, 9] and 
references therein, for studies in the Cartesian geometry), 
which constitute rotational flows. 

We find transformations between KP, cKP and ecKP 
equations, generalising the transformations between KP 
and cKP equations above, and use them to construct 
some important special classes of solutions of the derived 
version of the Kadomtsev-Petviashvili equation for both 
cases (i.e. for ecKP-I and ecKP-II). The ecKP equation 
is an integrable equation containing an arbitrary parame- 
ter a, and it reduces to the cKP equation both when this 
parameter tends to zero, and when r ^> a. The latter 
is the mathematical manifestation of the intuitively clear 
physical phenomenon: wave fronts will become nearly- 
concentric far away from the wave sources or boundaries, 
while they will "remember" their geometrical shape dur- 
ing the intermediate evolution. 

The KP and cKP equations belong to the family of 
universal integrable models of modern nonlinear physics 
(see, for example, [3 [10]). As well as for surface waves, 
these equations have been derived in many physical set- 
tings (see for example references in [5]), including in- 
ternal waves in a stratified fluid [TTJ [12] and, most re- 
cently, matter waves in Bose-Einstein condensates (BEC) 
(see [13] [14] and references therein) and cosmic dust-ion- 
acoustic waves [15]. Thus, the derived version of the 
equation could find many useful applications to the de- 
scription of the wave motion in problems where sources, 
boundaries and obstacles have elliptic or nearly-elliptic 
geometry. 

Our paper is organised as follows. In Section II we 
describe the derivation of the ecKP equation in the con- 
text of the classical surface gravity waves problem for 
an incompressible fluid. From mathematical perspective, 
the derivation for water wave problems is more challeng- 



ing than similar derivations for problems where nonlin- 
ear and dispersive terms are present in the equations, 
rather than originating from the free surface boundary 
conditions. In Section III we describe transformations 
between arbitrary solutions of the derived version of the 
KP equation and the original KP equation. Section IV is 
devoted to the lumps, line solitons and quasiperiodic so- 
lutions of the ecKP-I and ecKP-II equations. In Section 
V we discuss the approximate solutions for surface waves 
described by the derived equation. We conclude in Sec- 
tion VI by discussing the Lax pair of the derived equation 
and outlining possible applications and generalisations of 
our results. 



II. DERIVATION OF THE ELLIPTIC 
CYLINDRICAL KP EQUATION 

We consider the classical water wave problem for an 
incompressible fluid, described by the full set of Euler 
equations with free surface and rigid horizontal bottom 
boundary conditions . Since we aim to consider waves 
with the nearly-elliptic front, we write this set of equa- 
tions in the elliptic cylindrical coordinate system: 

x = d cosh a cos /?, 
y = dsinhasin /?, 

z = z, 

where the dimensional parameter d has the meaning of 
half the distance between the foci of the coordinate lines, 
and change the two horizontal components of the velocity 
vector appropriately: 

u — >• u cos j3 — v sin /?, 
v — > u sin /3 + v cos j3. 

Here, we keep the same notations (u, v, w) for the projec- 
tions of the velocity vector on the new coordinate lines. 

Let z = ho be the unperturbed fluid depth, A be 
the characteristic wavelength, p a the atmospheric pres- 
sure, and H the characteristic free surface elevation. We 
nondimensionalise the variables 



x — > Ax, y — > \y, z — > Hqz, t 




h^h + h 8 rj, p^p a + pg(h - z) + pgh p, 

which leads to the appearance of two usual nondimen- 
sional parameters in the problem: the long wavelength 
parameter 5 = and the small amplitude parame- 
ter e = ^j-, as well as a new nondimensional parame- 
ter 7 = which is n °t necessarily small. Scaling the 
dependent variables 



u — > eu, 



ev, w — >• ew, p 



ep, 



3 



we bring the full set of Euler equations in the elliptic 
cylindrical coordinates to the form: 



u t + e 



wu z + 



E{uu a + vup — v 2 ) + F(vu a — uup + uv) 



2je a G 



vt + e 



2je a G ' 

E(uV a + VVp + UV) + F(w a — UV/3 — u 2 



WV Z + 



E P(3 + Fp cx 
2-ye a G ' 



w t + e < ww z + 



2^e a G 



E(uw a + + F(vw a — uwp) 



2^e a G 

Z 1 

E(u a +vp+u) + F(v a - up +v) 
2^e a G 

P\z=l+eri( a ,f3,t) = T] -W e S 2 {f] aa + 77/3/3 



W z + 



(4) 
)' 

(5) 

(6) 
(7) 



G 



+ 



r] a sinh 2a + rjp sin 2/3 



2G 2 



i n I , , 2 <5 2 V* + Vp 



I 

3/2 ' 



(8) 



, ^ E(uri a +vrii3)+F(vri a -urii3) , , 
w\ z =i +ert ( a ,p,t) =Vt + e ^ , (yj 



^U=o o. 



(10) 



Here, W e = —^2 is the Weber number, and we denoted 

E = e 2a -cos 2/3, F = sin 2/3, 
G = sinh a + sin 2 /3. 

This set of equations reduces to the Euler equations 
written in cylindrical coordinates in the limit 

a — » 00, 7—^0 with 2^ eCK ~~ ^ r being finite. (11) 

The equation for linear waves (in the long-wave ap- 
proximation) is easily obtained from equations Q - ( [Tq] ) 
with e = 5 = as 



(12) 



7 2 (sinh 2 a + sin 2 /?) 
Note that equation (fl2|) indeed reduces to the equation 

(13) 



1 1 X „ 



for the long linear waves in the polar cylindrical coordi- 
nates in the limit ( pTj ). The derivation of the cylindrical 
KP (cKP) equation(also known as the nearly-concentric 
KP equation [2 , 4 ) is based on the existence of solutions 



of (|13|), which do not depend on /3, i.e. there exists an 
•f the equation ([IF 

Vtt ~ (Vrr + -Vr) = 0. 



exact reduction of the equation ([13]) to the equation 

1 



Unlike (13), equation (12) does not have an exact reduc- 



tion to the equation with no dependence on /?, which 
would seem necessary in order to derive a version of the 
KP equation for waves with nearly-elliptic front. Never- 
theless, such an equation exists as an asymptotic reduc- 
tion, and it turns out that this allows for a generalisation 
of the cKP equation to be derived. 
Next, we introduce the variables 



c 



5 2 



(7 cosh a — t) , 



R = —7 cosh a, 
o 4 



6 ■ R 

-~ snip, 



S 2 



U, v 



<5 3 



V, w 



<5 4 



V 



p 



tp 



which generalise a change of variables for the cylindrical 
coordinates [2]. Here, 27 cosh a is the nondimensional 
sum of the distances from a point on an ellipse to its 
foci. Thus, ( is an asymptotic characteristic coordinate 
for waves with nearly-elliptic front, and it becomes the 
characteristic coordinate for the concentric waves in the 



limit (hip. 

In tnese variables, the problem formulation Q - (10) 
assumes the following form: 



U c + P c + A [UU C + WU Z + Pr 
y/R 2 



R 



R 2 



-V c + 



A 2 
1 



A 2 

— (yP v + Rv 2 P c ) 



+ O(A 2 ) = 0, (14) 



VR' 2 - A 2 



+A 



uv c + wv z + 



R 



■ A 2 



A 2 



VR 2 - A 2 
(R-VR 2 -A 2 )(R 2 + A 2 

2(R 2 - A 2 ) 3 / 2 
(R-VR 2 -A 2 ) 2 + R 2 2 



vPr 



v'Pc 



2{R 2 - A 2 ) 
P z - /\W C + 0(A 2 



3/2 



U c + W z + A 
1 



U R 



VP„ 

= 0, 

R-VW 



+ (9(A 2 ) = 0, 



(15) 
(16) 



(Rv 2 U c + vUv) 



VR 2 - A 2 
+0(A 2 ) = 0, 



R? - A 2 
(R - VR 2 - A 2 )vV c + V„ + U 



P\ z =i + ah(c,r,») =H-AW e H cc + 0(A 2 



W\ z=1 +AH(C,R,r) 

W\ z =o = 0, 



-H c + AUH C + 0(A 2 



(17) 
(18) 
(19) 
(20) 
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where 

e 4 e 6 

Here, we have not shown the explicit form of the higher- 
order terms in the small parameter A (denoted by 
0(A 2 )) since these terms are not needed in the derivation 
of our asymptotic equation. 

We now seek an asymptotic solution of the system of 
equations and boundary conditions ( 14) - (20 ) in the form 



H = H + AH 1 +O(A 2 ), 

with similar expansions for J7, V, W and P. At leading 
order (0(1)) we obtain 



U 0C = Pi 



Vr 



1 



oc 



F-^Of H — vt ^0Ci 



VR 2 - A 2 VR 2 - A 2 

Poz = 0, U 0( + W 0z = 0, 
fb|*=i = Ho, WoU=i = "Hoc, ^o| z =o = 0, 

which yields, imposing the condition that the perturba- 
tion in U is caused only by the passing wave, 

P = H , U = H , W = -H 0( z, (21) 
1 „ R - VR 2 - A 2 



Vr 



oc 



VR 2 - A 2 



VR 2 - A 2 



-vH 0C .(22) 



At the next order (O(A)) we obtain the following equa- 
tions and boundary conditions: 



f/i C - Pk = p 0R + u u 0( + W U 0z 
R 2 _ A 2 \ vPo » + Rv P °<) 



(23) 



Vic 



VR 2 - A 2 



Plv ; — vP\C 

VR 2 - A 2 ( 



U V 0C + WoV 0z + R ,^ R2 . A -vPon 



(R-^R 2 -A*){R 2 + A 2 , _ 3 
2(R 2 - A 2 ) 3 / 2 



sjR 2 - A 2 
v 3 P 0( 



(R - VR 2 - A 2 ) 2 + R 2 2 

2(7? 2 _ ^2)3/2 V 

Plz = Woo 



(24) 
(25) 



U 1( + Wi 



F? _ A //?2 _ A 2 

-Uor + AO {Rv 2 U 0i + vU 0v ) 



R? - A 2 

L= ((7? - ^R? - A 2 )vV 0C , + V „ + U ) , (26) 



VR 2 - A 2 
PiU=i + H P 0z \ z=1 =H 1 -W e H 0(< , 
Wi\ z=1 + H W 0z \ z =i = -H K + U H 0< 
Wi|*=o = 0. 



(27) 
(28) 
(29) 



Then, <|2T}, j25) and (27) yield 

^2 



Pi 



-Hi 



0CC 



1 



■Hi, 



hence from (23), (26) and (29) we find, using (21), that 



iJi^ + 2Hqr + HqHq^ — — (2vHq v 

(R- VR 2 - A 2 )isV 0C 



+ 2Ris z H 0C ) + 
+VW + #o)U 



1 



VR 2 - A 2 



(30) 



Finally, substituting W\ into the remaining boundary 
condition (28), differentiating with respect to ( and us- 
ing ( [22] ), we obtain the elliptic cylindrical KP equation 
(ecKP) 



7/2 



C 



# 2 - A 2 



i? 2 - A 2 

How = 0. (31) 



Note that equation (24), written for the completeness of 



the set of equations, allows one to find V\ and is not used 
in the derivation of the ecKP equation. 
The scaling transformation 



with 



2a 



1/3 



7 = (Oa^a 2 ) 1 / 2 , 
(5 4^, 

a 2 = sign (a/3) 



brings the derived equation (31) to the form 
' Hr+mH( +H m + ^— rj H 
2 -H C ) + -^H„ 



12a 2 ci 2 (r 2 -a 2 Y 



C 



shown in the Introduction. Here, a = aA. If we let 
o; = l, then a — A and a 2 = sign Q — W e ). For typical 
water waves, a 2 = 1 (W e < |). However, a 2 = — 1 if 
the effects of surface tension are strong (W e > \)- It is 
natural to call the corresponding equations ecKP-II and 
ecKP-I, respectively, similarly to the terminology used in 
the Cartesian geometry. 
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III. TRANSFORMATIONS BETWEEN KP, CKP 
AND ECKP EQUATIONS. 

Considerations used to find the mapping from the solu- 
tions of the KdV equation to the class of solutions of the 
cKP equation [2] |4] can be extended to obtain transfor- 
mations between arbitrary solutions of all three versions 
of the KP equation, related to the Cartesian, cylindri- 
cal and elliptic cylindrical coordinates, respectively. The 
resulting transformations generalise the transformations 
between the KP and cKP equations [2j [6] , discussed in 
the Introduction. 

Indeed, the geometry of a wave with nearly-elliptic 
front, considered simultaneously in the Cartesian and el- 
liptic cylindrical coordinates, suggests the introduction 
of the sum and the difference of the nondimensional dis- 
tances from a point on the wave front to the two foci of 
the coordinate system 

d\ + d>2 = 27 cosh a, 
d\ — (I2 = 27 cos /?, 

where the foci have the following Cartesian coordinates: 
Fi(— 7, 0) and ^2(7, 0). We recall that the variables have 
been nondimensionalised, as discussed in Section II, and 
7 = |. Note that \{d\ + ^2) — t corresponds, up to 
the scaling, to the asymptotic characteristic variable 
introduced in Section II. 

Then, for the area satisfying ^3^, —> 0, we obtain 
the following asymptotic behaviour 

-(di+cfc) -t 

= \ (>/(* + 7) 2 + V 2 + y/{x - 7 ) 2 + y 2 ) - 1 

~x-t+-y 2 —— + 

4 \x + 7 x — 7 

Next, for sufficiently large a and small /3, our nondi- 
mensional variable x = 7 cosh a cos /? ~ ^-i?, and the 
previous asymptotics can be rewritten as 

I (dl + d a )-t~e+^^# 15> 

where ^ = x — t,Y = j^y and A = 7^ . Similarly, 
1 AY 2 

This asymptotic behaviour of the geometrically meaning- 
ful objects motivates the change of variables 



H (R,C,u)= V (R,Z,Y), 



where 



C = £+k 2 R 



Y 



2 i? 2 - A 2 ' 



V# 2 - A 2 



It is then verified by direct calculation that this transfor- 
mation maps the ecKP equation (31 ) to the KP equation, 
written in the form 



2f] R + 3t?% 



W, 



t]yy 0. 



To finish this section, let us summarise the transfor- 
mations between all three versions of the KP equation in 
water wave problems. We write the KP equation in the 
canonical form 



(U T + WU<: + U m ), + 3a 2 U YY = 0, 



(32) 



the cKP equation in the similar form 

Wt + q WWx + Wxxx + 1 w 

and the ecKP equation as 

' H T + 6HH ( + H ccc 



Za 2 



Wvv = 0, 



a 2 v 2 



2(t 2 - a 2 ) 

3<7 2 



12a 2 (t 2 -a 2 ) 
Then, the map 

W(r, X ,V) :=U[t, X 



H 

H uv = 0. (33) 



tV 2 
12a 2 ' 



-V 



transforms any solution of the KP equation into a solu- 
tion of the cKP equation, and the map 



U(t,£,Y)^ 
H(t,(,v) :=U[t,C 



12a 2 



(34) 
(35) 



transforms any solution of the KP equation into a solu- 
tion of the ecKP equation. Note, that the second trans- 
formation reduces to the first in the limit a — » 0. The 



map (35) also shows that for small a and small values 



of r any solution of the ecKP equation approaches some 
F-independent solution of the KP equation. These trans- 
formations can be inverted, and they can also be used to 
obtain the direct transformations between the cKP and 
ecKP equations. 



Indeed, the map inverting (35) has the form 
U(t,£,Y) :=h(t^ 



tY 2 



Y 



12a 2 t 2 



It transforms any solution of the ecKP equation into a so- 
lution of the KP equation. In particular, this map shows 
that for large values of r and finite values of Y any solu- 
tion of the ecKP equation will approach a y-independent 
solution of the KP equation (possibly, a constant or zero) . 



However, this limit is of academic interest only, since such 
large values of r are likely to lie outside of the range of 
applicability of the derived model, and we do not discuss 
this limit any more. 
The map 



W(t, X ,V)^ 
H(t,C,v) :=W[t,C 



2, ,2 



a v 



transforms any solution of the cKP equation into a solu- 
tion of the ecKP equation, and the map 



W(t, X ,V) ~H[t, X 



a 2 rV 2 



12a 2 (r 2 -a 2 )' 



V 



transforms any solution of the ecKP equation into a so- 
lution of the cKP equation. 







FIG. 1: (Color online) Solution to the ecKP-I equation ob- 
tained as the imag e o f the lump ( |36| ) with k = 1 under the 
action of the map (35) for a = 0.01 and several values of r. 



IV. SPECIAL SOLUTIONS OF ECKP-I AND 
ECKP-II EQUATIONS 

In this section we will consider some special solutions 
to the ecKP equation (33) per se, to illustrate the char- 
acteristic features of the equation. The considered ex- 
amples are exact solutions to the KP-I and KP-II equa- 
tions: lumps, line solitons and quasi-periodic solutions 
(see [T7J HE]), which become solutions to ecKP-I and 



ecKP-II under the mapping (35). 

If the a in (35) vanishes, the ecKP solution reduces 
to the corresponding cKP solution. For small "times" r, 
the cKP solutions look like solutions to the KdV equation 
(essentially no dependence on the transversal variable), 
whereas they develop horseshoe- type profiles for larger r. 
The ecKP solutions on the other hand show such profiles 
already for small r — a if a > 0. For large r (r ^> a) the 
solutions tend asymptotically to the corresponding cKP 
solutions. We will illustrate this behaviour at several 
examples. 

The first example we consider is the KP-I lump solu- 
tion, 



U{£,Y,t) 



4k(1 - k(£ - 3kt) 2 + k 2 Y 2 ) 
(1 + k(£ - 3kt) 2 + k 2 Y 2 ) 2 ' 



(36) 



with k = 1, under the mapping (35). It is visible that 
when t — y a and a is close to zero, the solution H ((", z/, r) 
is essentially independent of the coordinate v. This can 
be seen for a = 0.01 in Fig. [I] For larger values of a, the 
solution has a parabolic shape for small r — a as can be 
seen in Fig. [2] 

Next, we consider the 2-soliton solution of the KP-II 
equation in the form 

[/(£, Y, r) = 2dl In W (e^ 1 + , + ) , (37) 







FIG. 2: (Color online) Solution to the ecKP- I equation ob- 
tained as the image of the lump (36) with k = 1 under the 



action of the map (35) for a = 1 and several values of r. 



where 6j = kj^ + k^jY — 4k :|r, kj are arbitrary constants, 
and W is the Wronskian of the two functions. It can be 
seen for a = 0.01, i.e., close to the cKP case, in Fig. [3] 
where the formation of horseshoe waves can be clearly 
recognised. The corresponding ecKP solution for a = 1 
is shown in Fig. [4] where the curved profiles are already 
present for small r — a. 

Quasi-periodic solutions to the KP equations can be 
given in terms of multi-dimensional theta functions on 
compact Riemann surfaces of arbitrary genus n (see [19] ) 
in the form 

[/(£, F, r) = 2d 2 x In 9 [fr + Y v + rq + I] + C, (38) 
where is the Riemann theta function, p, v, q, I are peri- 
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FIG. 3: (Color online) 2-soliton solution ([37]) of the ecKP-II 
equation with for a — 0.01 with k\ — 1.5, — 0.5, k% — —2, 
k± = for several values of r. 



FIG. 5: (Color online) Genus 2 solution (|38) to the ecKP-II 
equation for a — 0.01 generated by the curve w 2 = Yli=i( z ~ 
ei), ei = -3, e 2 = -2,e 3 — -1, e 4 = 0, e 5 = l,e6 = 2 for 
several values of r. 




FIG. 4: (Color online) 2-soliton solution ([37} of the ecKP-II 
equation for a = 1 with k\ — 1.5, — 0.5, — —2, /C4 = 
for several values of r. 



ods of certain integrals on this surface, and C is constant 
with respect to the coordinates £, Y and r, see p~8j [21] 
for details. For a given Riemann surface and a given 
point on it, these quantities are uniquely determined. In 
genus 2 all such surfaces are hyperelliptic. In this case 
we consider the hyperelliptic curve with branch points 
— 1,-2,— 3, 0, 1,2. These solutions are numerically eval- 
uated with the spectral code by Frauendiener and Klein 
[20| |2T] . The related solutions to the ecKP-II equation 
are generated from the corresponding solutions of the 
KP-II equation via the mapping (35) with a = 0.01. We 
clearly see in Fig. [5] the formation of intersecting families 
of parabolic fronts. In the same setting with a = 1, i.e., 
a theta-functional solution to the ecKP-II equation, the 
formation of curved profiles is already present for small 



r — a as can be seen in Fig. [6j Both cases asymptotically 
coincide for r — > 00. 




FIG. 6: (Color online) Genus 2 solution ([38} to the ecKP-II 
equation for a — 1 generated by the curve w = Yl^=i ( z ~ e 0> 
e\ — —3, 62 — —2, C3 — — 1, €4 = 0, e.5 — 1, cq — 2 for several 
values of r. The formation of intersecting families of parabolic 
fronts can be clearly seen. 



In higher genus, the solutions are g-phase solutions, 
i.e., they have more structure as can be seen in Fig.[7J We 
consider here again hyperelliptic surfaces. The close to 
cKP solutions are for small time essentially independent 
of the transversal coordinate. 



FIG. 7: (Color online) Genus 3 solution (38) to the ecKP-II equation for a — 0.01 (left) and a = 1 (right), generated by the 
curve w 2 — Yl^i( z ~ e 0> e i — — 5, C2 — —4, e?, — —3, — —2, es — —1, e§ = 0, e7 = 1, es = 2 at r = a + 0.01. 



V. APPROXIMATE SOLUTIONS FOR 
SURFACE WAVES 

It is very interesting to see which approximate solu- 
tions for surface waves are described by the derived equa- 
tion. In order to do that we return to the original nondi- 
mensional variables x,y,t and re-parametrise our solu- 
tion as follows 

x = 7 cosh a cos /3, y = 7sinhasin/3, 



6 l/3 



(1 - 3W e ) 



1/3 



#(t,C,i/), 



where 



12a 2 (r 2 - a 2 ) 



2(r 2 -a 2 ) 
3a 2 



H 

H vv = 0, (39) 



C 



but rather numerical and real-life experimentation, and 
this is outside of the scope of our present paper. 

The 1-soliton solution of the ecKP-II equation is ex- 
plicitly written in the form 



K 2 



sech 



K 
2~ 



c 



TV 

~12 



-(K 2 +3L 2 )t + 5 )} 



(40) 



where K,L,5q are arbitrary constants. In the examples 
shown below we let the Weber number W e = 0, and 
the phase shift #o = 0- The corresponding surface wave 
elevation r] is plotted below for 7 = l,a = 2, A = 1/2. 
Similar solutions exist for the ecKP-I equation, as the 
image of an (unstable) line-soliton of the KP-I equation. 

The wave obtained when K = 1, L = is compact and 
symmetric, it is shown in Fig. [8] below. 



<j 2 = sign (1 — 3W e ) and r = R = a cosh a, 



c = 



7 A(l-3^)V3 ^ CQsh "- t )- 



Here, t is the physical time (nondimensional). Below, we 
assume that t > and consider the initial stages of the 
evolution. We also have e = a/^A and S = 2 A 3 / 2 . Since 

V a a 

asymptotic long-wave models usually provide a good 
qualitative (often quantitative) description even outside 
of the range of their formal asymptotic validity (i.e. the 
physical applicability of such models is usually wider than 
their formal asymptotic validity), we plot the solutions for 
all < j3 < 2tt and a > 0. Unless it is explicitly stated 
otherwise, at least parts of the shown solutions belong to 
the range of the formal validity of the asymptotic model 
(defined by r - 0(1), C - 0(1), 1/ ~ 0(1) as A ^ 0). 
The range of applicability of the model can not be de- 
fined within the scope of the formal asymptotic analysis, 







FIG. 8: (Color o nlin e) Surface wave corresponding to the one- 
soliton solution ( |4Q| ) of the ecKP-II equation with K — 1, L = 
for t = (top left), t = 0.25 (top right), t = 0.5 (bottom 
left), t = 1 (bottom right). 

For L / the solution has no symmetry with respect 
to the ^/-coordinate, as one can see in Fig. [9] for K = 1 
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and L = 0.1. The change L — »• — L yields the reflection 
of the wave front with respect to the x-axis 

rj(-L, x, y) = Tj(L, x, -y). (41) 

Therefore, it suffices to consider L > or L < 0. 








FIG. 9: (Color o nlin e) Surface wave corresponding to the one- 
soliton solution (40) of the ecKP-II equation with K — 1, L = 

0.25 (top right), t 



0.1 for t = (top left), t = 
left), £ = 1 (bottom right) 



0.5 (bottom 



When \L\ increases further, the compact nearly - el- 
liptic wave shape disappears. The wave becomes non- 
compact, and it rather describes the deformed line soli- 
ton, featuring an elliptic inhomogeneity in the central 
part of the wave. The solution is shown in Fig. [10] for 
K = l.L = -0.5. 







FIG. 10: (Color online) Surface wave corresponding to the 
one-soliton solution ( |40| of the ecKP-II equation with K — 
1,L = -0.5 for t = (top left), t = 0.25 (top right), t = 0.5 
(bottom left), t — 2 (bottom right). 

For sufficiently large \L\ solution is localised in the 
vicinity of some point satisfying the relation L sin fj = 
|L|, and strongly attenuates with time (Fig. 11). (Most 
likely, the large values of L lie outside of the range of 
applicability of the model.) 



FIG. 11: (Color onl ine) Surface wave corresponding to the 
one-soliton solution ( |40| of the ecKP-II equation with K — 
1, L = 2 for t = (leftJTt = 0.1 (right). 



In order to explain the observed features of the solution 



(40), let us note that the maximum of its amplitude is 
attained when the argument of sech 2 is equal to zero (if 
this is possible): 
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LVV 2 - a 2 v - (K 2 + 3L 2 )r = 0. 



This condition can be written either as 

2A 



2 - sin^ f3 



6 l/3 



[K 2 + 3L 2 ] cosh a 



+2A 1/2 6 1/3 Lsin/3sinha 



2t 



(42) 



or as 



sin 2 p - 2A 1/2 6 1/3 L tanh a sin j3 

+ -^r \ K<1 + 3^ 2 1 - 2 + — secha = 0. 

Let us first consider the case L = 0. For 

6 l/3 



(43) 



2A 



the solution has the form of a compact nearly-elliptic 
wave of narrowing width (as shown in Fig. [8|. Indeed, 

2A T - 2t , 2A o 

<- sech ^< 2 -^ K ' 

We note that although the width of the wave clearly 
changes, the amplitude is constant, which can be viewed 
as the manifestation of the solitonic nature of this solu- 
tion. 
For 



6 i/3 



6 1 / 3 1 
. <K 2 <— (^3.63 for A = -) 
2 A ~ ~ A v 2 ' 



the condition (43) is satisfied for 



sin/ 



and the wave splits into two deformed line solitons 
(shown in Fig. 12 ). 
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H 
n 1.2 





FIG. 12: (Color online) Surface wave corresponding to the 
one-soliton solution of the ecKP-II equation (40) with K — 
1.5, L = for t = (left), t = 2 (right). 



If 



K 2 > 



6 V3 



then the argument of (40) cannot be equal to zero, and 



the wave continuously attenuates with time. Most likely, 
these values of K lie outside of the range of applicability 
of the model. 

Similarly, we can consider the case L ^ (without loss 
of generality we assume that L > 0, see (41)). If 



6 l/3 



[K 2 +3L 2 ] >1, 



then the argument ( [42] ) of sech 2 is negative for L sin j3 < 
already at t = 0, and decreases with the increase of t. For 
L sin f3 = \L\ the argument is closer to zero than for other 
values of (3. Therefore rj has a maximum for L sin f3 = \L\ 
as a function of j3. Exact dependence of r] on a and /3 
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(or, equivalently, on x and y) can be seen in Fig 
However, as we mentioned above, such large values of K 
and L are likely to be outside of the range of applicability 
of the model, we mention them here for the completeness 
of our analysis, but we do not consider them any more. 
For sufficiently small K and L equation (43) yields 



sin/3 = A 1/2 6 1/3 Ltanha ± VD, 



where 



D = 2 — - A6 1/3 L 2 sech 2 a - — secha. 
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Then, for a — >• oo we obtain 



sin (3 = A 1 / 2 6 1/3 L± \ 2 



2A 

6V3 



K 2 . 



Therefore there exist several natural domains for the 
values of the parameters K and L. Let us note that 

2A 9 J I 

if 2 1/3 K < 0> tnen tne som tion (|40|) continuously 



attenuates, and we do not consider this range of values. 
2A 9 

Let 2 trK a > 0. Then, there are several cases. 



(i) If 



A l/2 6 1/3 L± 2 



2A 

^173 



K 2 



> 1, the solution (40) 



has the form of a deformed nearly-elliptic wave (shown in 



Fig. [9]). Indeed, in this case the argument of (40) can not 
be equal to zero for sufficiently large a, and therefore the 
solution is localised. For sufficiently small L the above 
inequality implies: 



A X / 2 6 1/3 L- 

A x / 2 6 1/3 i- 
which yields 



2A 



K 2 > 1 and 



2A 

6V3 



K 2 < -1, 



2 -^ >1 + A l/2 6 l/3 Li 



and 



K + 3L 2 



g2/3 gV3 

AV2 < ^A"- 



(44) 



It is easy to see that if (44) holds, then for any /3 
2A 



2 - sin 2 p - ^ [K 2 + 3L 2 ] > 2A 1 / 2 6 1 / 3 L \sm(3\ . 



gl/3 

Therefore we can introduce the notations 
2A 

2A l/2 6 l/3 isin( g = A ^ sinh ( ao ) ) 



2 - sin 2 (3 



[K 2 + 3L 2 ] = A((3) cosh(ao), 



where 



A(f3) = { 2 -sin 2 /? 



2A 

6~T73 



[X 2 + 3i 2 



^AG^^sin 2 /?} 1 / 2 , 

2A 1 / 2 6 1 / 3 isin ( 9 
tanhcto = 



2 - sin 2 $ 



2A 

6V3 



[K 2 + 3L 2 } 



In these notations equation (42) takes the form 
cosh(o; + ao) = 



A(/3) 7 



(45) 



Therefore in this case for any /? and for sufficiently large 
t there exists the value of the parameter a such that the 
argument (42) will be equal to zero. Because all func- 



tions are continuous and different iable, the maximum 
of H(t, C, v) will be attained along some smooth closed 
curve. For sufficiently small t and for any a, f3 the value 
of the function H(t,(,v) will be less than K 2 /2. Since 
the sign of ao depends on the sign of /?, for L ^ it 
follows that the solution will be asymmetric with respect 
to (3. For Lsin/3 < the wave will be wider than for 

L sin (3 > (shown in Fi g. [9|. 

/ OA 

(ii) If A 1 / 2 6 1 / 3 L + J 2 - —K 2 > 1 and 1 > 



6 i/3- 



^/ 2 &/*L- 



2A 

6V3 



i^ 2 > — 1, the nearly-elliptic wave 
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breaks for negative y, producing a single deformed line- 
soliton (shown in Fig 10). 



(ii) if 



A V2 6 1/3 L± 2 



2A 

6V3 



K 2 



< 1, the nearly- 



elliptic wave breaks both for negative and positive y, pro- 
ducing a pair of deformed line-solitons (shown in Fig. 13). 





FIG. 13: (Color online) Surface wave corres pon ding to the 
one-soliton solution of the ecKP-II equation ( 40 1 with K = 
1.6, L = 0.1 for t = (left), t = 2 (right). 

Apart from these generic cases, there are also some 
exceptional cases (corresponding to the crossings from 
one generic case to another). To illustrate that, we show 
a solution corresponding to the condition 



2A 

6V3 



K 2 = 1, 



(46) 



two 



which is the borderline case in between the last 
generic cases. The solution is shown in Fig. 14 , and it 



can be interpreted as describing the splitting of the wave 
looking like a resonant Y-soliton (or 'Miles soliton' [16]), 
known from the theory of the original KP equation, into 
two deformed line solitons. The existence of such solu- 
tions might indicate the instability of the Y-soliton with 
respect to perturbations in the area of the crossing. 





FIG. 14: (Color online) Surface wave corresponding to the 
exceptional one-soliton solution of t he ecKP-II equation ( 40 ) 
with K = 1.5 and L defined by |46f (L « 0.1) for t = (left) 
and t = 1 (right). 

Let us note that the one-soliton solution of the KP-II 



equation used to obtain the solution (40) of the ecKP-II 



equation above can be written in the form of the 'canon- 
ical' soliton for the KP-II equation 



U{t,Z,Y) 
ipj = exp(/cj^ 



2<9f ln(^i 



-^2), where 



■ak)Y 



4*£r) 



if K = ki — &2 and L = (k\ + k^)a (we let Sq = 0). Then, 
the conditions on kj follow from the conditions on K 



and L, discussed above. We show in Fig. [15] the solution 
of ecKP-II equation which corresponds to the canonical 
one-soliton solution of KP-II equation with a = l,fci = 
0.5,^2 = -0.4. 





FIG. 15: (Color online) Surface wave corresponding to the 
solution of the ecKP-II equation obtained from the canonical 
KP-soliton with ki = 0.5, k 2 = -0.4 for t = 0.5 (left), t = 1 
(right). 

Using Darboux transformations, one can obtain the 
canonical two-soliton solution of the KP-II equation in 
the form 

U(t, e, Y) = 2d\ ln(0i0 2 £ - 2 0ie), where 

(up to the phase shifts, which can be added to the 
phases). The corresponding two-soliton solutions of the 
ecKP-II equation are shown in Fig. 16 and Fig 17 (asym- 



metric and symmetric two-soliton nearly-elliptic waves, 
respectively) . 







FIG. 16: (Color online) Surface wave corresponding to the 
two-soliton solution of the ecKP-II equation with k\ — 
0.5, = -0.4, k 3 = 0.4, U = -0.3 for t = 1 (top left), 
t = 2 (top right), t = 3 (bottom left), t = 4 (bottom right). 

Finally, let us choose W e = 2/3 and consider the ecKP- 
I lump solution 



if(r,0) 



4k(1 - k(C + tv 2 /\2 - 3kt) 2 + k 2 (t 2 - a 2 )v 2 ) 



(1 + k(C + tv 2 /12 - Zkt) 2 + k 2 (t 2 - a 2 )v 2 ) 2 
It is easy to see that for sin j3 = the wave elevation r\ 
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FIG. 17: (Color online) Surface wave corresponding to the 
two-soliton solution of the ecKP-II equation with k\ = 
0.5, k 2 = -0.5, k 3 = 0.4, fc 4 = -0.4 for t = (top left), 
t = 0.5 (top right), t = 1 (bottom left), t = 2 (bottom right). 




FIG. 18: (Color online) Surface wave corresponding to the 
ecKP-I lump solution with n — 0.25, for t — (top left), 
t = 0.5 (top right), t = 1 (bottom left), t = 2 (bottom right). 



has the form 



-A- 3 Jlf( C -3KT), 



where the function 



f(X) 



4«(1 - kX 2 ) 
(1 + kX 2 ) 2 



has one high maximum for X = and two weak minima 
for X = zLa/3/k, where 

/(0)=4k, f(±y/3fi) = -k/2. 

Therefore, for sufficiently large values of £ the wave 
elevation 77 has two deep minima 

16ft fa 
Vmin = -Qijsy- 



for 



sin /3 = 0, cosh a = 



The corresponding surface wave elevation rj is plotted 
in Fig [18] for 7 = l,a = 2, A = 1/2, W e = 2/3 and 
ft = 0.25 . 



VI. CONCLUDING REMARKS 



The ecKP equation ( 39 ) derived in our paper is an in- 



tegrable model, which can be obtained as a compatibility 



condition of the following linear problem (Lax pair): 

aVv — V r2 — a2 ^CC + (V 7-2 — a 2 H(r, C, v) 

j- 22 
tQ a v 



12Vr 2 - a 2 144a 2 Vr 2 - a 2 
-4^ccC" (6IT(t,C,i/)- 



(47) 



12a 2 (r 2 - a 2 ) 



3^ c (r,C,^) + 



3aif(r,C,^) 
Vr 2 - a 2 



a 2 C^ 



12a(r 2 - a 2 ) 3 / 2 



^ (48) 



Indeed, the compatibility conditions have the form 



2(r 2 - a 2 
2 



H 



12a 2 (r 2 - a 2 ) ' t* - a 

When a = we recover the Lax pair of the cKP equation 
(e.g. 0). 



A solution ^ ec xp(r, z/) of the linear system ( |47| ), (48) 
is expressed via the solution ^xp(t, £,y) of the linear 
system of the KP equation (U^ = Uy) 

^ r = -4^-6C/(r,e,^ 
-(3C^(r,e,r) + 3(7£r(r,e,y))^ 



as follows 



ipecKp(rX,v) = ^kp r, C 



x exp < 



Qtv 



TV r- 

^ 3 (r 2 + a 2 ) 



12a Vr 



432a Vr 2 - a 2 



Let us also note that the functions H and J7 are related 

by 



TV 



H(t, C, i/) = Vt 2 - a 2 U (r, £, F) - ^C/ (r, £, F) . 
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The Jacobian of the transformation 

f = C - tv 2 I(YI(j 2 ), Y = VV 2 - a 2 */, 

used to construct some special solutions of the ecKP in 
the preceding sections, is equal to \fr 2 — a 2 . It is positive 
and bounded for all r > a. (This condition is automati- 
cally satisfied for our derivation.) 

In our paper we considered only some simple solutions 
of the derived equation. Recently, there has been signif- 
icant progress in the classification of soliton solutions of 
the KP equation with applications to water wave prob- 
lems (see [22H25] and references therein). It would be 
interesting to see the counterpart of this classification for 
the derived equation, and for the approximate solutions 
for surface waves. 

The derivation of the ecKP equation from the full set 
of Euler equations opens the way to the study of internal 
and surface waves on a current for a stratified fluid, as 
well as accounting for the effects of variable background 
and Earth's rotation, which will be reported elsewhere. 
It paves the way for other applications, for example, in 
the context of matter waves in Bose-Einstein condensates 
(e.g., P~3j[l4]), since the hydrodynamic form of the Gross- 



Pitaevskii equation is similar to the problem formula- 
tion (TJ. Also, recent studies of 'spherical nebulons' [26] , 
based on the spherical KP equation, can be extended 
since the ellipsoidal KP equation, associated with the el- 
lipsoidal coordinates, can be derived from the equations 
for 'a dusty plasma' along the lines discussed in this pa- 
per. 

Finally, it is natural to ask a question whether one can 
derive other versions of the KP equation, associated with 
other coordinate systems (i.e. with other wave geome- 
tries), and whether one can find the general description 
of all admissible maps of the type discussed in Section 
III, associated with the problem formulation ([!]). 
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